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Abstract 

The construction of the Zhang-Rice singlet is revisited in the light 
of recent understanding of high-temperature superconductors at optimal 
doping. A minimal local model is derived which contains the physical 
regime found relevant for ARPES experiments, characterized by signifi- 
cant direct oxygen-oxygen hopping. For the values of orbital parameters 
indicated by experiment, the Zhang-Rice singlet is strongly mixed with 
a pure oxygen singlet of the same symmetry. The destabilization of the 
Zhang-Rice ground state is due to the oxygen singlet having twice as 
large a coherence factor with respect to oxygen-oxygen hopping. An anal- 
ogous quantum phase transition is identified in the t-t'-J model. The 
orbital-antisymmetric copper-oxygen singlet is confirmed to be irrelevant, 
as found originally. The usual symmetry analysis is extended to include 
dynamical symmetries. 

PACS: 74.72-h 74.25. Jb 71.10.Hf 

1 Introduction 

The nature of the electronic wave functions which undergo the transition to su- 
perconductivity in copper-oxide perovskites is still unknown. It is generally at 
least tacitly assumed that the superconducting (SC) transition at optimal dop- 
ing does not invoke radically different states than those already present at the 
Fermi level slightly above T c . In this approach, shared here, the understanding 
of the normal state is considered the key to the SC mechanism. 

The antiferromagnetic (AF) instability in the vicinity of half-filling in the 
same compounds is by contrast well understood [I] as a charge-transfer ef- 
fect, dominating the first (large Mott-Hubbard interaction Ud) Coulomb cor- 
rection [5] to the ideal (paramagnetic) Mott insulator [3] . Theoretical attempts 
to connect this underdoped region with the SC region on the hole-doped side 
abound in the literature. They have received early support in the ionic limit 
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by the Zhang-Rice [I] construction, which considered an added hole in an AF 
background, when double occupation of the copper site was blocked by a large 
Ud- It is known that the hole then occupies the oxygen sites, each bridging two 
coppers in the CuC>2 plane. Assuming that the hole prefers to bind with one of 
the coppers, Zhang and Rice have shown that a uniform spreading of the hole 
around that copper, in an orbitally symmetric "Zhang-Rice singlet" (ZRS) is 
energetically strongly preferred, and propagates like a projected fermion. This 
provided a way to vindicate the single-band t-J model as a physical reduc- 
tion of the Emery (three-band Hubbard) model [SJ [3] , at least for underdoped 
systems, as an alternative to the usual formal derivation, which applies the 
Foldy-Wouthuysen transformation to the one-band Hubbard model [7J. 

Challenges to this picture come from two directions. One is fundamental, 
that there is no a priori reason for the hole to prefer one copper over the other. 
When it is undecided, triplet states are strongly admixed, while it may still 
move about as a fermion [§J. The other is to accept the basic idea of a single 
Cu04 molecule, but consider additional states and overlaps within it, to find 
the limits of the particular construction employed [5]. In the present work, the 
natural next step in the Zhang- Rice construction is taken, which introduces just 
one additional singlet, occupying the oxygen sites. 

As will be shown below, such an extension is the minimal model to ex- 
press the low-energy physics of optimally doped high-T c superconductors in 
the small-cluster language. At optimal doping, the Zhang-Rice ground state be- 
comes strongly admixed with said oxygen singlet, chosen by the same symmetry 
arguments as originally used by Zhang and Rice [4]. 

The simplicity of the calculation is used as an opportunity to introduce 
classification of state vectors by dynamical symmetry on a small example. This 
method may have some potential for larger problems. 

2 Physical scales 

Quantum chemical calculations on small clusters predict the regime 

U d > A pd > \t pd \ > \tpp\, (1) 

used by Zhang and Rice, and confirmed in XPS experiments [12] . Here A pd = 
£ P — £d is the copper-oxygen splitting, positive in the hole picture, while t pd 
and t pp are the copper-oxygen and oxygen-oxygen overlaps, respectively. In the 
limit HI), the ZRS is the ground state. The energy scale of the XPS experi- 
ments and corresponding calculations is 5-10 eV. In high-T c materials, there 
is a universally acknowledged crossover scale of the order of 0.1 eV between 
the lowest and higher energies, although it is still debated whether its origin 
is phononic [13] or magnetic [13]. Neither small-cluster nor LDA calculations 
obtain this crossover, and the latter are consistent with the high-energy result 
(Q}. A recent analysis [15] showed three distinct many-body scales in BSCCO 
which were not reproduced by contemporary LDA calculations, the lowest of 
which was the above-mentioned crossover, the middle the so-called "waterfall," 
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not of interest here, and the highest an "anomalous enhancement" of the LDA- 
predicted band-width, by about 0.8 eV visible below ep in optimally doped 
BSCCO, i.e. roughly twice that in full. In principle, LDA calculations are able 
to reach lower energy scales than small clusters. Thus no ab initio effort has so 
far reached the vicinity of the Fermi level in high-T c superconductors. Below 
the crossover scale, the physical regime of the electrons actually undergoing the 
transition to superconductivity must still be inferred from phenomenological 
analysis of low-energy experiments. 

The optimally doped region is equally well connected to the metallic over- 
doped region as to the underdoped one. Attempts to understand the high-T c 
perovskites at optimal doping starting from the paramagnetic metal have a long 
standing on the hole-doped side [10] . They have mostly been framed in terms of 
the mean-field slave-boson (MFSB) formalism, allowing for effective multiband 
approaches, in which the values of the orbital parameters depart significantly 
from those expected on the basis of high-energy experiments and LDA calcula- 
tions [TT] . 

Instead of taking the high-energy "bare" regime as the starting point, and 
introducing many-body effects to reach lower energies, MFSB phenomenology 
goes in the opposite direction. Beginning with unprejudiced zeroth-order fits at 
the Fermi level, one introduces the minimal perturbation necessary to model the 
lowest crossover. Such an analysis of angle-resloved photoemission (ARPES) [TH] 
has recently been extended to electron-doped compounds [T7] , and produced the 
interesting result, that their lowest-energy sector is in a very similar physical 
regime as in the hole-doped ones [16j HB] : a large direct oxygen-oxygen hopping 
t pp , significantly larger than the effective copper-oxygen hopping t p d, and with a 
band- width of the same order as the effective copper-oxygen splitting A p d > 0, 
in the hole picture: 



The ARPES scales on which this observation is based are 0.2 eV in BSCCO and 
0.5 eV in NCCO. The lowest-energy crossover was modelled by a {n,ir) boson, 
whose scales were associated with AF in the same compounds. The low-energy 
regime ([2]), markedly different from ([!]), is characterized in particular by an 
"anticrossing" of the wide oxygen band with the copper level, so that the Fermi 
level is found in an oxygen-dominated, strongly dispersive part of the Brillouin 
zone, even though Ed < e p . As will now be shown, a directly analogous regime 
exists in the ionic limit. 

3 Model and technique 

The Hamiltonian of a single CUO4 molecule is 



U d > i\t pp \ >e p -e d = A pd > t 2 pd /A pd . 



(2) 
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Figure 1: Orbital parameters and symmetries used in the article. In this sign 
convention t pp < is physical [T§] . 
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where P a = p_ rj . a + p_~ g — p +£ a — p + ~ a describes the four oxygens around 
the copper, and hda, h p .i<j are the number operators of the copper and oxygen 
orbitals. The relevant orbitals and parameters [TI5] are shown in Figure Q] All 
25 states of two holes with opposite spins on the five orbitals are taken into 
account. 

Standard symmetry analysis is used to reduce the Hamiltonian to block- 
diagonal form. The basic symmetry is of the square, corresponding to the 
group -D4. Take the z-axis, orthogonal to the molecular plane, to be the spin 
quantization axis. Since the spins have total projection zero, flipping all spins 
is a symmetry operation, and since all orbitals are planar, this is equivalent to a 
reflection in the plane. Then D^h is the magnetic (color) group of the problem, 
and both the Zhang- Rice singlet, \ZR), and the doubly-occupied copper state, 
|d 8 ), fall in its A\ u representation, containing five singlets. 

As usual, the main symmetry group does not fix the state vectors in full. It 
is easy to complete the reduction by hand, but instructive to identify the reasons 
for the remaining ambiguities. One is, as expected, symmetry with respect to a 
subgroup, Di h c D±h- The particular D2/1 group involved contains the spin- flip 
and three 180° rotations, about the z-axis and about the two diagonals at ±45° 
relative to the Cu-0 axes. It stands out because its operations change the sign 
of the copper-oxygen overlap t p d, while the energies depend only on t^ d . 

The more interesting reason is dynamical symmetry, imposed by the as- 
sumed degeneracy of the four oxygen orbitals. It eventually affects only two 
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Bi u vectors, but to fix them "honestly," the full dynamical symmetry appara- 
tus must be invoked, even on such a small problem as the present one. The 
group elements are the Hamiltonian terms which break the oxygen degeneracy. 
Two are visible in ©, proportional to +t pp and — t pp , respectively. If these are 
represented by matrices in one-particle space, multiplying them generates the 
third, which corresponds physically to next-nearest-neighbor (n-n-n) hops from 
one oxygen to the other on the same axis. These three (and the identity) form 
the Abelian group D2, appearing here as a dynamical symmetry group. The 
n-n-n term must be diagonalized simultaneously with the sum of t pp terms in 
([3]) to separate out the last vector without manual intervention. 

The outstanding physical feature of dynamical symmetries is that they can 
connect states with doubly and singly occupied degenerate sites, obviously im- 
possible by geometrical transformations. Two technical advantages of dynamical 
over geometrical symmetry groups are that they grow with the size of the prob- 
lem, and depend on topology rather than shape. Both are significant when the 
clusters are motivated by solid-state investigations, and contain many equiva- 
lent sites. E.g., the dynamical group of the four oxygens in two-particle space 
is Z|, with sixteen elements, again Abelian. It is not needed here because the 
geometrical symmetry already reduces the problem to very small blocks. On 
the other hand, as a molecule becomes larger, the blocks due to dynamical 
symmetry are expected to become smaller than those selected by geometrical 
symmetry, especially if the molecule is of odd shape, as is common in Monte 
Carlo calculations. 



4 Result 

The complete reduction gives ten stationary states unaffected by t pd . In addition 
the Hamiltonian contains six 2x2 blocks and one 3x3 block. The latter is found 
in the A\ u subspace, connecting the states |cZ 8 ), \ZR), and an oxygen two-hole 
singlet, \p~ 2 ): 

\P- 2 ) = J(|l) + |2> + |3) + |4) (4) 
+ | 12) - |13) - |14) - |23) - 1 24) + |34)) . 

where \i) = are doubly occupied states, and \ij) = — |j|) 

are antisymmetrized singly occupied states. The numbers 1-4 refer to the four 
oxygens, as in Fig. [T] The block reads: 

\d s ) : / 2e d + U d 2V2t pd \ 

\ZR) : 2V2t pd e d + £ p + 2t pp 2^2t pd (5) 
\P~ 2 ) ■ \ 2V2t pd 2e p + 4t pp J 

The states \p~ 2 ) and |c? 8 ) have an equally large coherence factor 2\/2t pd with 
respect to \ZR), as found by perturbation theory [4]. However, the coherence 
factor of the oxygen-oxygen overlap t pp is twice as large for \p" 2 ) as for \ZR). 
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Figure 2: Parameter space of the model (J5J). 



Since the physical sign of tpp is negative, this is a generic mechanism to desta- 
bilize the ZRS for all values of Ud (including Ud — > oo), which identifies \p~ 2 ) 
as a relevant perturbation of the Zhang-Rice ground state. 

By the same analysis, the orbital-antisymmetric singlet |0^4), also considered 
by Zhang and Rice, appears in a 2x2, Bi u block with another oxygen singlet, 
call it \pp): 

\OA) : f e d + e p -2t pp 2t pd \ , . 

\pp) ■ \ 2t pd 2e p ) [0) 

The splitting here is also A p d + 2t pp , as between \p~ 2 ) and \ZR) in ([5]), but 
it occurs because the copper-oxygen singlet is antibonding (t pp < 0), while the 
oxygen singlet is non-bonding with respect to t pp . Its structure is 

\PP) = ^= (-11) + 12) ~ |3) + |4> + |13) - |24)) . (7) 

The orbital-antisymmetric state remains irrelevant in the presence of oxygen- 
oxygen hopping. The tp^-coherence factors in ([5]) and ([6]) may be identified 
by perturbation theory [1], without the need to distinguish between \p~ 2 ) and 
\pp). However, once t pp is introduced, the different effects on \ZR) and \OA) 
are related to the different symmetry properties of \p~ 2 ) and \pp). The Zhang- 
Rice singlet \ZR) is destabilized because \p~ 2 ) has a large coherence factor 
4£j, p , while \OA) moves upwards "out of the way," and \pp) cannot influence the 
ground state because it is unaffected by t pp (but of course it can anticross with 
\OA) as the latter rises). 

The coherence factors in the reduced model (O are the largest in the whole 
space, so for \t p d\ > and t pp < it contains the ground state. The relevant 
parameters are Ud/ 'A pt ; and t pp / A p d- The precise value of t p d is irrelevant, as 
long as one is in the parameter regime ([2]), corresponding to optimal doping. 
Choosing it low for graphical reasons, one obtains Figure [21 At the origin 
(Ud = 0, t pp — 0) is the weak-coupling regime, where the ground state is 
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dominated by |d 8 ). As Ud increases, this state goes through two anticrossings 
with the others, ending up highest in the Zhang-Rice regime (Ud large, \t pp \ 
small), where \ZR) is the ground state. Now as \t pp \ increases, the state \p~ 2 ) 
goes down at twice the rate of \ZR), so it must eventually anticross with it. 
Thus when both Ud and \t pp \ are large, the ZRS is no longer the ground state. 

5 Discussion 

The present paper describes a simple and robust symmetry-based mechanism 
which prevents the isolation of the ZRS as the only candidate for the ground 
state of the CUO4 ion when Ud is large. To wit, if the oxygen-oxygen hopping t pp 
is significant, the local environment will be in a charge transfer regime involving 
two states, the ZRS and an oxygen singlet. The two are mixed by the symmetry- 
enhanced overlap 2^t p d across the symmetry-reduced gap A p d + 2t pp . This is 
the metallic counterpart to the result [I] that the antiferromagnetic interaction 
J near half-filling is also dominated by the charge-transfer term t^ d / A pd , while 
the Mott- Anderson contribution t^ d /U is negligible. 

It has long been established that ARPES fits require significant ratios of 
\tpp/t p d\ in order to agree with the data [TT]. Recently there appeared detailed 
ARPES fits, both qualitatively and quantitatively in accord with experiment, 
in practically the same orbital regime for NCCO [17] and BSCCO [16] , which 
established the anticrossing regime ([2|) as an appropriate phenomenology for 
these compounds, meaning the ratio \t pp \/A p d is also enhanced, especially so on 
the electron-doped side. 

Changes in orbital parameters as a function of doping are a standard feature 
of effective theories such as mean-field slave-boson (MFSB), which account for 
the effects of strong correlations by changing t pd and A p d- The question, as al- 
ways with effective band models, is what this means for the local environment: 
do local hops "see" the bare or the renormalized band parameters. The conclu- 
sion that the effective band regime at optimal doping is reflected in an analogous 
local regime, which implies the destabilization of the Zhang-Rice singlet by the 
above symmetry mechanism, is supported by two independent arguments, one 
phenomenological and one theoretical. 

The recently emphasized |15) failure of LDA calculations to reproduce exper- 
imentally observed band widths is phenomenologically significant, because LDA 
predicts a smaller width than is actually observed, so any conceivable inclusion 
of additional correlations would only make the discrepancy larger. Notably, 
LDA calculations agree with high-energy experiments and calculations in the 
ionic limit. The anomalous enhancement of hopping via the oxygen sites is con- 
sistent by symmetry with the observed discrepancy, which is large in the nodal 
direction, where the oxygen band is strongly dispersive. The enhancement of 
the ARPES bandwidth of 1-2 eV over LDA predictions [15 is quantitatively in 
good agreement with the value of t pp needed to enter the anticrossing regime, 
i.e. destabilize the ZRS. 

Theoretically, the parameter space of the three-band MFSB model has been 
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extensively investigated P~8] from the point of view of outcomes, namely: which 
regions of renormalized orbital parameters are indicative of strong renormaliza- 
tion of input parameters by the MFSB mechanism, and which imply preservation 
of the bare parameters. In the so-called "resonant-band" regime, characterized 
by a narrow effective bonding band and copper-dominated wave functions, the 
quasiparticle weight in the conduction band is significantly reduced, being pro- 
portional to the doping S instead of the total hole content 1 + S. In that case, 
the band calculations cannot be directly related to a local picture, since the 
doped-hole weight is locally either one or zero. However, in the anticrossing 
regime @, the opposite occurs: MFSB theories with Gaussian fluctuations de- 
part from the resonant-band regime and acquire a full quasiparticle weight in 
the Brillouin zone, only with a reduced t p d [18] • They are therefore comparable 
and consistent with the present ionic ones, as long as both are in the anticrossing 
limit. Thus the present work reaffirms the possibility [20] that the optimizations 
leading to a large direct oxygen-oxygen overlap below the crossover scale occur 
already at the local orbital level. 

Realistically, t p d is not so small, so the ground state in real materials retains 
substantial copper content. In the local picture, this means a significant ZRS 
content in the ground state. However, even if t p d is as large as \t pp \, giving 
a large anticrossing gap ~ A p d, the ground state still contains nearly equal 
components of \ZR) and |p~ 2 )- The system is physically in the charge-transfer 
regime, which also provides the dominant contribution lj to J. In that case, 
the reduction to one band by the limit [4] A p d = Ud — A p d — * oo is purely 
formal, since that limit replaces the charge-transfer regime a priori by the Mott 
regime. 

While the MFSB calculations allow both for regimes with strong renormal- 
ization and no renormalization [IS], the derivation of the t-J model from the 
ZRS [4] implies by construction that the band parameters are locally relevant. 
However, it has been noticed [9] that in order to face experiment, the t-J model 
needs to be supplemented by additional next- and next-next-nearest neighbor 
hops, making up the so-called t-t'-t"-J model. Since the t' term breaks particle- 
hole symmetry in the same way as the t pp term in the three-band model, its 
appearance in practical calculations [21] is immediate evidence that the oxygen 
degree of freedom cannot be eliminated from the physical picture near opti- 
mal doping. The need to introduce t' at the band level translates microscop- 
ically into the need to admix the oxygen singlet ([5]) to the ZRS, to describe 
the local electronic environment around optimal doping. (When the oxygen is 
treated explicitly in the band picture, t pp has the same sign for both n- and 
p-doping [16j [17] , as also found in LDA calculations [22] , consistently with its 
simple chemical interpretation.) 

To see the real effect of the t' term, the dispersion of the t and t' terms may 
be rewritten 

— 2i(cos k x + cos ky) — At' cos k x cos k y = 

-4(i + t') + 4(i + 2t')(sin 2 *E + sin 2 - 16t' sin 2 ^ sin 2 ^L, (8) 
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separating contributions by particle-hole symmetry. The second term is man- 
ifestly symmetric around half-filling, while the third is not, since it achieves 
its minimum value (zero) along the non-dispersive lines k x = and k y = Q. 
In practice tt' < 0, and the ratio \t'/t\ — 0.2-0.4. For t' = — 0.3i, a com- 
mon enough parametrization around optimal doping, |4£'/(f + 2t')| = 3, so 
the symmetry-breaking term enters with three times the contribution of the 
symmetry-preserving term. Thus the apparently small values of \t' /t\ found in 
the literature are quite misleading in terms of the actual amount of particle-hole 
symmetry breaking needed to achieve realistic Fermi surface fits within the t—t'— 
J model. In particular, for t' = —t/2 the symmetry-preserving term vanishes, 
so the t-t'-J model passes through a quantum phase transition in the vicinity 
of its applicable parameter range. 

By symmetry, this is the same type of transition appearing in the present 
paper as the destabilization of the ZRS in the three-band case. Physically, they 
are different insofar as the three-band model destabilizes the ZRS by oxygen 
hopping in low order, while in the t-t'-J model the Zhang-Rice singlets are 
preserved by assumption, and the transition occurs by the additional t' overlap 
among them, i.e. in high order in the original three-band model. But to account 
for ARPES data, both models agree that such a transition is approached near 
optimal doping. 

The limitation of the above arguments to the metallic phase should be kept 
in mind. Near half-filling, the oxygen hopping is effectively reduced by the anti- 
ferromagnetic interactions, and precisely at half-filling, the system is insulating 
independently of the value of t pp . In the insulating phase, the high-energy ionic 
limit is recovered, both theoretically and experimentally. As t pp becomes more 
effective with doping, locally the first admixture to the ZRS appears, identified 
here as the oxygen singlet ([5]). Thus a broad picture of the optimal-doping 
regimes of high-T c compounds emerges, in which NCCO is deep in the anti- 
crossing regime, i\t pp \ >• A p d, and BSCCO is at its limit, 4|£ pp | w A p d- On the 
other hand, the remarkable Fermi-surface evolution of LSCO with doping [33] 
is indicative of strong renormalizations of the band parameters, especially when 
the system is underdoped. In such a resonant-band regime, it was also shown 
directly [21] that the strong doping-dependence of t pp in the MFSB fits did not 
imply a similar variation of the local parameters. 

The present analysis points to one pitfall of considering the oxygen part of 
the ZRS as a single-particle state in its own right, as if the hole on the cop- 
per were unhybridized with the oxygens. However, the same caution applies 
to the oxygen singlet found here. In particular, one should be careful in inter- 
preting the corresponding Bloch states. The propagation of strongly correlated 
cluster states in a lattice is at present largely unexplored in more than one di- 
mension, with controlled analytic results available only in special limits, which 
themselves amount to uncontrolled simplifications. Such is the limit of infinite 
dimensions [35], physically corresponding to diffusion, or the Falicov-Kimball 
limit [26], which is dispersive, but replaces dynamical constraints by geometri- 
cal ones. 

The model presented here is as sparse as possible, allowing for the exper- 
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imentally observed regime with only three state vectors. To avoid being too 
minimal, some attention should be paid to the neglected Coulomb integrals U p , 
V pp and V v d-, the on-site oxygen, intracell oxygen, and copper-oxygen repulsions, 
respectively. They do not break the symmetries above, introducing only matrix 
elements within each symmetry block. In particular, they connect the ground- 
state block ([3]) with the other two A\ u states by matrix elements Up/A, smaller 
than the terms 2\f2t p d- On the diagonal of this block, the Coulomb terms can 
be reabsorbed in the site energies and copper on-site repulsion Ud- For realistic 
V pp <C Up, this leads to a small reduction in the copper-oxygen splitting A p d, 
and practically no change in Ud, while for large V pp = Up/ 2, a useful limiting 
case, the effective Ud is increased and the effective A p d decreased, both by about 
1 eV. Thus the inclusion of neglected integrals mildly tends to push the system 
further into the anticrossing, Ud — > oo regime. 

An extension of the Zhang-Rice argument to the electron-doped side is only 
possible with larger molecules, containing more than one copper atom. This 
would allow to study the competition of magnetic and charge ordering at the 
local level. It also opens a way to compare the Zhang-Rice and Emery- Reiter [5] 
constructions directly on the hole-doped side. Obviously, the spaces involved 
are much larger. It is on these larger problems, both theoretical and numerical, 
that hope is invested in the assignment of quantum numbers by dynamical 
symmetries, as introduced in the present work. 

6 Conclusion 

It has been found that in the presence of local repulsion on copper, the di- 
rect oxygen-oxygen overlap is a relevant perturbation of the Zhang-Rice singlet. 
With the possible exception of LSCO, real hole- and electron-doped high-T c 
superconductors appear to be in a physical regime where it is strongly mixed 
with a pure oxygen singlet. 
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